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Abstract: Future space-based experiments, such as OWL [1] and JEM-EUSO [2], view large atmospheric and terres- 
trial neutrino targets. With energy thresholds slightly above 10^® eV for observing airshowers via air fluorescence, the 
potential for observing the cosmogenic neutrino flux associated with the GZK effect is limited. However, the forward 
Cherenkov signal associated with the airshower can be observed at much lower energies. A simulation was developed to 
determine the Cherenkov signal strength and spatial extent at low-Earth orbit for upward-moving airshowers. A model 
of tau neutrino interactions in the Earth was employed to determine the event rate of interactions that yielded a tau lepton 
which would induce an upward-moving airshower observable by a space-based instrument. The effect of neutrino attenu- 
ation by the Earth forces the viewing of the Earth’s limb to observe the rv-induced Cherenkov airshower signal at above 
the OWL Cherenkov energy threshold of ~10’^®’® eV for limb-viewed events. Furthermore, the neutrino attenuation lim- 
its the effective terrestrial neutrino target area to ~3 x at eV, for an orbit of 1000 km and an instrumental 

full Field-of-View of 45°. This translates into an observable cosmogenic neutrino event rate of ~l/year based upon two 
different models of the cosmogenic neutrino flux, assuming neutrino oscillations and a 10% duty cycle for observation. 
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1 Introduction 

Future space-based air fluorescence experiments employ 
wide field-of-view optics from a orbiting platform(s) to 
monitor a vast amount of the atmosphere. For the OWL 
(Orbiting Wide-angle Light Collectors) mission, the mass 
of the viewed atmosphere corresponds to more than 10^^ 
metric tons (mtons). The design choices for OWL were 
driven by the goal to measure the UHECR spectrum, via 
the air fluorescence technique, with high statistics above 
^10^® eV. Studies indicated that the ability to measure neu- 
trino interactions in the atmosphere via the air fluorescence 
signature exists, but the predicted event rate based upon 
cosmogenic neutrino flux models [3] is < 1/year (assum- 
ing a duty cycle of 10%) due to the paucity of neutrino flux 
above 10^® eV. Furthermore, the neutrino event rate quickly 
diminishes as the energy threshold becomes further away 
from the ~3 x 10^® eV threshold for full neutrino aperture, 
which assumes both OWL satellites stereoscopically view 
each event. 

Airshowers also produce an intense, beamed Cherenkov 
signal. OW’L simulation studies indicated that the energy 
threshold for observing the optical Cherenkov signal form 
a nadir- viewed, upward-moving vertical airshower initiated 
by a particle at sea level would be 10^^-^eV. OWL also 
views a large terrestrial area; assuming 1000 km orbits and 


the two OWL satellites are tilted to view a common area, 
the terrestrial area monitored ranges from 6 x 10° km^ for 
a 500 km satellite separation to nearly 2 x 10® km^ for a 
2000 km satellite separation. These vast areas offer a large 
neutrino target, depending upon the depth that provides a 
measurable signal. Tau neutrino charged-current (CC) in- 
teractions offer a mechanism to exploit this large, terrestrial 
neutrino target; at high energies, the produced tau lepton 
has a sufficient Lorentz-boosted length to escape the Earth, 
decay in the atmosphere, and create an upward moving air- 
shower that could be observed via the Cherenkov signal. 
This paper details the calculations used to quantify the sen- 
sitivity of measuring the cosmogenic neutrino flux using 
space-based measurements, assuming the performance de- 
fined by the OWL experiment, of the Cherenkov signal cre- 
ated from upward-moving tau neutrino induced airshowers 
originating in a terrestrial neutrino target. 

2 Optical Cherenkov Signal Simulation 

In order to determine the optical Cherenkov signal strength 
and profile at a orbiting instrument, a computer-based sim- 
ulation was constructed based upon parameterizations de- 
scribed by Hillas [4, 5]. The charged-particle density for 
an airshower is given by the Greisen parameterization as a 
function of shower age, with the airshower electron angu- 
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Figure 1 : Schematic of the Cherenkov signal generated by 
an upward-moving airshower induced by a tau neutrino in- 
teracting in the Earth. 

lar distribution and the energy-dependent, charged-particle 
track length fraction defined by well-behaved mathemati- 
cal functions dependent upon shower age. The index of 
refraction of air is given as a function of atmospheric gram- 
mage, or altitude, which defines the Cherenkov angle and 
Cherenkov energy threshold. The latter, combined with 
the charged-particle track length fraction, determines the 
fraction of the particles in the airshower that contribute to 
the Cherenkov signal for a given altitude. The atmosphere 
is described by the Shibata parameterization [6], and the 
wavelength-dependent attenuation of UV light in the atmo- 
sphere is described by a parameterization [7] based upon 
more detailed models. 

The parameterizations are then used to generate the 
Cherenkov signal at an arbitrary altitude for an upward- 
moving airshower. The lateral spread of the charged par- 
ticles in the shower was not considered since this is a 
small effect for large viewing distances. A 100 m air- 
shower step size was used starting at sea level, and the sam- 
pling of the electron angular distribution at each step was 
< 10“^ radians. The charged particle track length frac- 
tions were sampled in decades of energy for each step, from 
the Cherenkov energy threshold to a factor of 10 below the 
total airshower energy. The Cherenkov light was gener- 
ated at each step from 200 nm to 600 nm in increments 
of 25 nm, but the effects of ozone absorption, in the at- 
mospheric attenuation model, significantly reduces the sig- 
nal below ~300 nm. Two numerical azimuthal integrations 
were performed: one defined by the Cherenkov angle about 
the electron angular sampling vector, and the other by an- 
gular sampling vector about the shower direction. 

The effects of the Earth’s curvature were modeled as these 
become important for viewing angles away from the nadir 
direction. Figure 1 illustrates the geometry of this effect in 
relation to an upward moving airshower induced by a tau 
neutrino interaction in the Earth near the Earth’s limb. The 
Cherenkov airshower simulation accounts for this effect by 
determined the proper angle, shown as ip in the figure, at 
each shower propagation distance step, relative to the view- 



Figure 2; The radial profile of the simulated Cherenkov 
light density at an altitude of 1000 km for upward- moving 
10^^ eV airshowers viewed at angles of 0°, 45°, and 55°. 

ing angle Oy ■ Thus the atmospheric attenuation depth and 
the distance from the shower and the UV detector in orbit 
can be accurately determined. 

The results of the Cherenkov airshower simulation are 
shown in Figure 2 for Oy = 0°,45°, and 55° at an alti- 
tude of 1000 km. The curves show the Cherenkov light 
density (photons/m^) as a function of radial distance from 
the projected location of the airshower core in the plane 
perpendicular to the airshower direction at 1000 km alti- 
tude. The form of the curves show that the Cherenkov light 
cone (illustrated in Figure 1), is approximately uniform in 
photon density up to a radius which closely corresponds to 
that defined by the Cherenkov angle at shower maximum 
and the distance from shower maximum to the measure- 
ment. This implies that a good energy resolution can be 
achieved by sampling the uniform part of the distribution, 
e.g. by one of the two OWL satellites. The lateral size 
of the Cherenkov light cone also defines the solid angle: 
1.3 X 10~^ steradians for the nadir {6y = 0°) case and 
1.9 X 10“^ sr for 9y = 55°, which corresponds to view- 
ing close to the Earth’s limb (at an altitude of 1000 km, the 
Earth’s limb is given by 9y w 60°). A comparison of the 
results presented in Figure 2, after scaling the energy and 
altitude, are in good agreement (30% difference) with re- 
sults which employ a full airshower Monte Carlo and more 
detailed atmo.spheric attenuation modeling [8]. 

The factor of ~10 decrease in the photon density for Oy = 
55°, as compared to the nadir case, is mainly due to the 
effect of the Earth’s curvature. A factor of ~5 reduction 
is due to a 1/r^ effect (for 0° the distance from 1000 km 
altitude to shower max is 990 km while it is 2190 km for 
the 55° case) and a factor of '^2 reduction is due to atmo- 
spheric attenuation of the Cherenkov signal. 

The signal strength in an OWL ’eye’ is defined by the 3 
m optical aperture, the optical transmission, and the quan- 
tum efficiency of the focal plane detector. When these 
are combined, the photo-electron sienal streneth in an 




32nd International Cosmic Ray Conference, Beijing 201 1 


Er 

(eV) 

Depth (km) 

(P = I7S) 

Depth (km) 
(P = 2.3^) 

Mass (mtons) 
for 10® km® 

10^° 

0.05 

0.05 

1.3 X 10^4 

IQie 

0.5 

0.5 

1.3 X 10^® 


5 

5 

1.3 X IQi® 

IQis 

29 

16 

7.6 X 10^® 


18 

10 

4.8 X lOi® 


Table T. The terrestrial neutrino target depth for Vr CC 
interactions, for water and rock targets, and the target mass 
for 10® km^ for rock as a function of energy. 


OWL pixel is approximately the photon density in units 
of photons/m^, which unit of the ordinate axis in Figure 2. 
Thus a Cherenkov signal of 400 photons/m^ corresponds 
to a signal in a OWL instrument of ^400 photo-electrons. 

3 Cosmogenic Neutrino Rate Determination 

For charged-current tau neutrino interactions in the Earth, 
the target depth is estimated by the propagation distance of 
the tau lepton: 7 ct. At energies above 10^® eV, neutrinos 
and antineutrinos have similar cross-sections for neutrino- 
quark scattering, and the produced lepton carries ~75% of 
the incident neutrino energy at 10^® eV with the energy 
fraction increasing to ~80% at 10^® eV [9]. Er = is 
assumed for this ICRC paper. At higher energies, catas- 
trophic energy losses (a convenient parameterization is 
used [10]) limit the effective depth of the neutrino target 
above ~10^* eV. Table 1 details the terrestrial neutrino tar- 
get depth for interactions, as a function of E„, for both 
water and rock and calculates the mass for a target area of 
10® km^ assuming rock {p = 2.3 g/cm®). Note that tau en- 
ergy losses limit the target depth above ~10^® eV while the 
maximum neutrino target mass also occurs at ~10^® eV. 
Assuming a terrestrial neutrino target area of 10® km®, the 
rate of cosmogenic neutrino interactions that lead to an 
observable tau-induced, upward-moving airshower via the 
Cherenkov signal can be calculated by numerically inte- 
grating the neutrino flux, cross-section, target mass, . . ., 
with respect to the neutrino energy, and accounting for the 
solid angle of the Cherenkov signal. Two cosmogenic flux 
models were used, the Bartol (Oa = 0.7) model [3] and 
a more recent calculation by Scully&Stecker [11]. The 
models, both which assume proton dominated UHECRs, 
provide the electron and muon neutrino and antineutrino 
fluxes. These were summed for each model and divided by 
3 to predict the Vr flux, eg assuming neutrino oscillation 
yields equal flavor fractions at Earth. 

Figure 3 shows the results presented as observable 
events/year per energy decade for the two cosmogenic 
models, assuming a solid angle of 1.5 x 10”® sr, a terres- 
trial neutrino target area of 10® km® and a 100% duty cycle 
for the experimental measurement of the Cherenkov signal. 
The rate prediction based upon the Bartol flux is given by 
the upper curve while that based upon the Scully&Stecker 
model is given by the lower curve. An interesting feature is 
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Figure 3: The rate of observable tau neutrino interactions 
via the airshower Cherenkov signal for two cosmogenic 
neutrino flux models. The rate is presented in number of 
events/year per energy decade. 

that the event rate does not become appreciable until above 
10^® eV, indicating that the neutrino target mass, and, to 
a lesser extent, the neutrino cross-sections, becomes suffi- 
ciently large at this energy to yield a significant rate. 

While the neutrino rate prediction of > 100 events/year is 
encouraging, the calculation does not take into account the 
fact that the duty cycle for observation of the Cherenkov 
signal is closer to 10% and the Earth is relatively opaque to 
these energetic neutrinos. This latter point forces the orbit- 
ing instrument to view the limb of the Earth to observe air- 
showers induced by neutrinos with sufficiently short path 
lengths in the terrestrial target. Figure 4 details the geome- 
try: the orbiting detector must be tilted to view the limb of 
the Earth, and the effective terrestrial neutrino target is con- 
strained to be no larger than an appropriately short neutrino 
path length in the Earth. 

The problem now effectively becomes one of analytic ge- 
ometry. The tilt angle of the instrument is defined by its 
EOV and the angle to the Earth’s limb, the calculation of 
the viewed elliptical area on the Earth’s surface is deter- 
mined by the intersection of the tilted cone defined by the 
FOV and a sphere of radius 6378 km. The 2-dimensional 
(eg flat) area of the ellipse can be determined using the 
conic section relationship between the angle of the plane, 


Altitude (km) 

Full Area (km®) 

Effective Area (km®) 

500 

1.7 X 10® 

1.6 X 10® 

1000 

3.8 X 10® 

2.6 X 10® 

2000 

8.9 X 10® 

3.9 X 10® 


Table 2: The full and effective viewed surface areas for an 
45° FOV instrument, tilted to view the Earth’s limb, as a 
function of altitude. The effective area is constrained by 
chord equal to the v interaction length in Earth at 10^^ eV. 
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Figure 4: A 2-dimensional schematic of the effective neu- 
trino target area constrained by the p interaction length. 
Only those neutrinos with a chord length sufficiently small 
will be relatively unattenuated by the Earth. While a tilted 
wide FOV UV imager monitors a large, elliptical area of 
the Earth, only a small (darkest region) portion samples a 
significant neutrino flux. 

which defines the ellipse in the cone, and the eccentricity of 
the ellipse [12]. The 2-dimensional area of the viewed el- 
lipse was then increased by 5% to account for the effects of 
the sphere’s (Earth’s) curvature. The effective area (dark- 
est area in Figure 4) is defined by the chord whose length is 
the neutrino attenuation length in Earth, for a particular en- 
ergy. This describes a truncated ellipse which is inscribed 
in a rectangle of calculable dimensions. Assuming the trun- 
cated ellipse can be approximated by a parabola of similar 
width, one can use the fact that the area of the parabola is 
2/3 that of the bounding rectangle [13] to approximate the 
area of the truncated ellipse, which is the effective area of 
the terrestrial neutrino target. 

Table 2 presents the total, tilted viewed terrestrial area for 
an instrument with 45° FOV and the effective neutrino tar- 
get area at 10^^ eV, as a function of altitude, determined 
from the analytic geometrical calculation. While the full 
area monitored by a tilted instrument is substantial, the ef- 
fective neutrino target area available near the Earth’s limb 
is reduced by more than a factor of 10, to 2.6 x 10^ km^, 
assuming a 1000 km altitude. 

Combining the energy-dependent effective area results 
within the numerical integration using the two cosmo- 
genic flux models and assuming 10% duty cycle, the pre- 
dicted observable tau neutrino event rate is reduced to ~2 
event/year using the Bartol model and ~1 event/year using 
that of Scully&Stecker. 

4 Conclusions 

Simulation studies of the Cherenkov signal from upward- 
moving airshowers indicate that an orbiting experiment 
with a 3 m optical aperture and UV sensitivity of an OWL 


instrument would have an energy threshold of slightly 
higher than 10^® eV for airshowers generated near the 
Earth’s limb. This is well-matched for tau-induced air- 
showers generated by cosmogenic neutrino interactions in 
the Earth. However, the 10% duty cycle inherent for the 
Cherenkov observation and the Earth’s attenuation of the 
neutrino flux limit the effective terrestrial area, which is 
estimated to be ~3 x 10® km^ for E^, — 10^^ eV, as- 
suming an instrument with a 45° full FOV tilted to ob- 
serve the Earth’s limb. Using two different cosmogenic 
neutrino flux models, the predicted observable event rate 
is ~1 event/year. While factor of 2 improvements may 
be available using different orientations of the two OWL 
satellites or realizing a gain in the duty cycle, the net ef- 
fect of these improvements may be balanced by potential 
decreases caused by more realistic modeling of the energy 
distribution of the create tau lepton, the airshower gener- 
ated by the tau decay, and the inherent shower fluctuations. 
While the tau-induced airshowers could be observed via 
the air fluorescence technique, which has a much larger ob- 
servational solid angle than that inherent to the Cherenkov 
signal, the higher energy threshold of 10^® eV (and a 
factor of ~10 higher for viewing near the Earth’s limb), 
severely limit the sensitivity to the cosmogenic neutrino 
flux. This reinforces a similar result from more detailed 
Monte Carlo studies of the ability of OWL to measure air- 
showers induced by cosmogenic electron neutrinos in the 
atmosphere, which predict < 1 event/year. However, if the 
energy threshold for air fluorescence could be reduced to 
10^® eV for an OWL-type mission, studies have indicated 
that the measurable cosmogenic neutrino event rate would 
be ~50 events/year for the atmospheric channel, assum- 
ing a 10% duty cycle. There could also be a significant ob- 
servable rate from interacting in the Earth observed via 
air fluorescence, if the energy threshold could be reduced. 
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Surface Area ~ 10® km^ 
SatSep=500 km: 5.8 x 10® km^ 
SatSep=2000 km: 1.6 x 10® km^ 




Motivation 
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a' with 45° full FOV optics, 3 m optical aperture 
els with ps readout, 0.06° pixel angular size 
irbits, stereo configuration 
icence Energy Threshold ~ lO^^ eV 
=^ 10 % 
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Upward Airshower Cherekov Model 




Airshowers generate a strong Cherenkov signal 

Model of Upward-going (relative to Earth) created based upon 
parameterizations given by Hillas (JPhysG 8) 

• Index of refraction of Air vs altitude 

• Track length fraction vs shower age and energy 

• Definition of E^her/ sample in decades for E > E(~her 

• Angular distribution of particles vs angle, energy, and age 

Use to build Cherenkov light generation model based upon Greissen 
shower parameterization in 100 m steps (200 nm < X < 600 nm| 

• Model atmospheric attenuation 

• Shibata atmospheric model (see Gaisser's book) 

• Rayleigh scattering (see Sokolsky's book) 

• Ozone absorption (Krizmanic 28^^ ICRC, Salt Lake City) 

• Effects of Earth's curvature modeled 


• Aerosols not modeled (est of signal reduction at large angles == 
25% for 1.2 km scale height and X^=14 km ) 
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100 PeV Airshower Cherekov Results 
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Vertical Airshower 



Nadir Pointed Instrument 

Vertical Airshower (0v=O°) 
Signal: 400 photons/m^ 

PE Signal: = 400 PEs in OWL pixel 
Lateral Size: ~ 40 km at 1000 km 
Solid Angle: 1.3 x 10'^ steradians 
Distance from 990 km 


Tilted Pointed Instrument 

Slanted Airshower (0^=55°) 


PE Signal: ~ 35 PEs in OWL pixel 
Lateral Size: = 110 km at 1000 km 
Solid Angle: 1.9 x lO'^ steradians 
Distance from Sj^g^: 2187 km 
55° at z=1000 km -> 71 deg at z=0 




-30 -20 -10 0 10 20 30 

Distance from Shower Axis (km) 


-100 -80 -SO -40 -20 0 20 40 60 80 100 

Distance from Shower Axis (km) 


Slanted signal reduction: % factor due to larger distance (1/r^) 
and 60% additional reduction due to atmospheric scattering. 
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Consider Interactions: {Gandhi et al., PRD 58} 

• (vbarIM) = 5.5 x 10"^® (E^ / 1 GeV)®-^®^ ^ > ^gie 

® <E^> ~ 0.75 X at = 10^^ eV, rising to == 0.8 x at = lO^o eV 

Depth of neutrino target given by the minimum of: 

• yet for the produced taon 
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Cosmogenic Neutrino Flux and v Interaction Rates 


0^ESST 



Total Events/Year 

^BARTOl./'S = 204 
^Scully+Steckar/ 3=1 25 


Bartol: Engel et al., PRD 64 
Scully & Stecker; AstroPartPhys 34 



Vt- 'Event Rates' per year per decade in 
Assuming: 

1. 10^ km^ area 

2. target depth prescription 

3. Q = 1.5 X 10'^ steradians 




5. 100% duty 
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Geometry of the Calculation 


Instrument is rotated to view limb of Earth at edge of FOV 

• Viewed surface area is ellipse on a sphere 
Effective area limited by transparency of the Earth to v 

• Effective area is defined by the chord determined bv ' 


Iffeciive area is defined by the chc 
a ns pa re II c:y: - 1421 km at 10^^ 

id decreases as E increases 
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1.7 X 10« 

1.6 X 10® 

3.8 X 10® 

2.6 X 10® 

8.9 X 10® 

3.9 X 10® 
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Cosmogenic Event Rates vs Effective Area 


Including the energy-dependent Neutrino attenuation by the Earth 
the duty cycle leads to a more realistic Event Rate determination 


■■ . 



Relevant Parameters: 

® 10% Doty Cycle 

* 1000 km altitude 

® 45® full Field"Of“View 

* eV for less} Airs 
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Conclusions 

• Upward-moving airshowers are in principle observable from space-based 
experiments, with an energy threshold dependent upon the size of the optics: 

• for OWL parameters, ~ 10 PeV, depending on viewing angle. 

• While a large area of the Earth is viewed by OWL-type instruments, the effective 
area for "Earth-skim mirig' neytrinos is limited fay the attenuation of the neutrino 

signal by the Earth 

• Geometric calculations demonstrate the observable rate for the cosmogenic 
neutrino flux via observation of the Cherenkov airshower signal is - 1 event/year, 

assuming 10% duty cycle. 

• However, these results re-inforce those obtained by studies of interactions in the 
atmosphere observed via the air fluorescence signal: if the energy threshold for the 
air fluorescence measurements can be significantly reduced (eg to 10^^ eV), the 
event rate predicted by cosmogenic flux models would be appreciable: 

* this would require at least a factor of 3 increase in the size of the optics: 
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